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PURPOSE:  The  Coastal  Engineering  Technical  Note  (CETN)  herein  describes  the  Numerical 
Model  of  the  Longshore  current  for  Current  and  Waves  (NMLong-CW)  that  accounts  for  the 
interaction  between  a  current  and  surface  waves.  NMLong-CW  can  simulate  wave 
transformation,  the  steady-state  wave-generated  longshore  current,  and  change  in  water  level  at 
an  inlet  by  waves  and  wind.  The  interaction  between  a  current  and  waves  can  significantly  alter 
the  wave  height  and  wavelength.  NMLong-CW  is  a  one-dimensional  model  and  is  limited  to 
situations  were  longshore  uniformity  applies. 

BACKGROUND:  The  original  numerical  model  NMLONG  (Kraus  and  Larson  1991;  Larson 
and  Kraus  1991)  was  developed  under  the  U.S.  Army  Corps  of  Engineers’  Dredging  Research 
Program.  It  calculates  nearshore  wave  transformation,  water  level  change,  and  longshore  current 
across  a  single  profile  line,  under  the  assumption  of  longshore  uniformity  in  the  profile  and 
hydrodynamic  processes.  NMLONG  solves  the  wave  energy  flux  conservation  equation, 
accounting  for  wave  shoaling,  refraction,  breaking,  and  reforming.  The  transformation  of 
random  waves  is  calculated  by  Monte-Carlo  simulation.  Wave  energy  dissipation  accompanying 
depth-limited  breaking  is  described  in  accordance  with  Dally,  Dean,  and  Dalrymple  (1985).  The 
cross-shore  momentum  equation  is  solved  numerically  to  obtain  the  change  in  water  level,  and 
the  equation  governing  alongshore  momentum  gives  the  distribution  of  the  steady  (time  average 
on  a  wave  scale)  longshore  current  velocity  across  the  shore.  The  longshore  current  and  change 
in  water  level  (setup  and  setdown)  driven  by  a  local  wind  can  also  be  simulated.  Nonlinear 
bottom  friction  is  computed  from  an  efficient  approximation,  and  lateral  mixing  is  modeled  with 
an  eddy  viscosity  that  depends  on  the  local  wave  orbital  velocity  and  wave  height. 

Both  tidal  and  wind- generated  currents  can  be  comparable  to  or  exceed  the  strength  of  the  wave¬ 
generated  longshore  current.  Currents  produced  independently  of  waves  as  by  wind  and  the  tide 
may  be  in  opposite  directions,  producing  complex  distributions  of  the  current  across  the  shore. 
The  capability  of  representing  the  action  of  currents  in  NMLONG  and  the  interaction  between 
the  current  and  waves  resulted  in  a  substantially  revised  model  and  name  as  NMLong-CW, 
where  CW  stands  for  interaction  between  currents  and  waves.  In  the  revised  model,  wave 
transformation  in  the  presence  of  a  current  is  calculated  from  conservation  of  wave  action  flux 
instead  of  wave  energy  flux,  and  the  wave  dispersion  relationship  is  modified  to  include  a 
current. 

NMLong-CW  was  revised  to  describe  both  depth-  and  steepness-limited  wave  breaking,  and  a 
unified  formulation  to  quantify  energy  dissipation  produced  by  these  two  types  of  breaking  is 
employed.  As  before,  the  longshore  current  and  mean  water  level  are  calculated  through  the 
alongshore  and  cross-shore  momentum  equations,  respectively.  However,  in  the  alongshore 
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momentum  equation,  provision  for  an  external  (large-scale)  current  has  been  incorporated  in 
addition  to  the  wave-  and  wind-driven  currents.  At  present,  the  external  current  must  be 
specified  as  an  input,  although  in  the  future  an  option  will  be  added  allowing  the  user  to  generate 
current  distributions  within  NMLong-CW  corresponding  to,  for  example,  large-scale  tidal 
currents  and  ebb  jet  flows.  Additional  background  information  on  the  current-and-wave 
interaction  can  be  found  in  other  CETNs  (Smith  1997,  1999)  produced  under  the  Coastal  Inlets 
Research  Program  (CIRP). 

In  the  following,  a  summary  is  given  of  the  equations  used  in  NMLong-CW  with  focus  on  the 
enhancements  made.  The  procedures  for  calculating  wave  transformation,  longshore  current, 
and  water  level  change  are  discussed  separately.  Also,  a  section  is  included  on  wave  blocking 
and  the  criterion  applied  to  describe  this  phenomenon.  Capabilities  of  NMLong-CW  are 
demonstrated  by  examples.  The  model  is  operated  through  a  graphical  interface  that  runs  on  the 
Windows  95/98  and  NT  platforms  for  personal  computers.  NMLong-CW  is  at  the  state  of  the  art 
in  calculation  of  nearshore  waves  and  currents,  and  this  CETN  documents  the  underlying  physics 
implemented. 

WAVE  TRANSFORMATION:  Under  the  condition  of  alongshore  uniformity,  wave 
transformation  across  a  nearshore  profile  is  described  by  the  equation  for  conservation  of  wave 
action  flux  (e.g.,  Jonsson  1990): 


d  (  ECga  cos  (3  N 
dx  ^  co;. 


Ed 

(0,. 


(1) 


where  E  is  the  wave  energy  (linear  theory  used),  Cga  the  absolute  wave  group  speed,  [5  the  wave 
ray  direction,  co,  the  relative  wave  frequency  (=  2n/Tr,  where  T,  is  the  relative  wave  period),  ED 
the  wave  energy  dissipation,  and  x  the  cross-shore  axis  pointing  offshore.  The  presence  of  a 
current  with  magnitude  U  and  direction  8  will  alter  the  wave  transformation,  and  the  expressions 
for  Cga  and  /j  may  be  obtained  from  geometric  considerations  to  yield  (see  Ligure  1  for  a 
definition  sketch): 
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where  Cgr  is  the  relative  group  speed,  and  a  the  direction  of  the  wave  orthogonal.  The 
definitions  of  the  angles  are  -90  <  a  <  90  and  -180  <  8  <  180. 
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Figure  1.  Definition  sketch  for  waves  propagating 
on  a  current  (overbar  implies  a  vector) 


To  determine  the  wave  properties  at  a  certain  total  depth  d=  h  +  /),  where  h  is  the  still-water 
depth,  and  77  is  the  mean  deviation  from  still-water  level  (set  up  or  set  down),  the  dispersion 
relationship  including  a  current  has  to  be  solved  (Jonsson,  Skougaard,  and  Wang  1970): 


J—tanhkd  = 

u 


j~d~  4  U  cos(S-a )Ta  d 

V4l "  d  L ) 


(4) 


in  which  Ta  is  the  absolute  wave  period,  L  the  wavelength  at  depth  d,  k  the  wave  number 
(27t/L),and  La  denotes  the  deepwater  wavelength  neglecting  the  current  (=gTa2/2n).  Equation  4  is 
solved  by  a  Newton-Raphson  technique,  which  gives  rapid  convergence.  The  relative  group 
speed  is  determined  from: 


C 


gr 


(1+-20-) 

^  sinh  2 kd  J 


(5) 


and  the  relative  phase  speed  C,  is  computed  as: 


C  = 


—  tanh  kd 
k 


(6) 


in  which  g  is  the  acceleration  of  gravity.  As  can  be  seen  from  Equation  4,  the  wave  angle  a  must 
be  known  before  the  wave  properties  can  be  calculated.  Thus,  Snell’s  law  is  employed  to 
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determine  wave  refraction  and  the  variation  in  a  across  the  profile.  For  a  wave  traveling 
between  locations  denoted  by  indices  1  and  2,  Snell’s  law  is  written: 

sincfj  _  sin  a,  ^ 

I  I  ’ 

Equations  4  and  7  must  be  solved  simultaneously  because  both  a  and  L  are  unknown  at  the  next 
grid  point  in  the  numerical  calculations. 

The  wave  energy  dissipation  in  Equation  1  is  given  by  (compare  Smith,  Resio,  and  Vincent 
1997): 

ED=^-(E-EJClr  (8) 

dD 

where  K  is  an  empirical  coefficient  (=  0.15),  dD  the  length  scale  controlling  the  dissipation  (equal 
to  the  water  depth  d  in  Dally,  Dean,  and  Dalrymple  (1985)),  and  Es  the  energy  of  a  stable  wave 
for  which  breaking  ceases  and  a  wave  can  reform.  To  generalize  Equation  8  to  describe  all  water 
depths  and  both  depth-  and  steepness-limited  breaking,  the  following  expressions  are  employed 
for  the  stable  wave  energy  and  the  dissipation  length  scale: 
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where  p  is  the  density  of  water,  T  is  an  empirical  coefficient  (=  0.4),  and  y b  is  the  ratio  between 
wave  height  and  water  depth  at  incipient  breaking  under  depth-limiting  conditions  (typically 
taken  to  be  0.78).  The  wave  height  at  incipient  breaking  is  calculated  from  the  Miche  criterion, 
modified  by  Battjes  and  Janssen  (1978)  to  be  applicable  for  all  water  depths: 
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In  shallow  water,  Equation  12  reduces  to  Hi,  =  y ),d,  indicating  that  Equations  8-11  recover  the 
original  formulation  by  Dally,  Dean,  and  Dalrymple  (1985).  Equation  12  is  applied  with  the 
local  water  depth  to  determine  Hi,  from  which  Hs  is  obtained  from  Equation  10. 

WAVE  BLOCKING:  Waves  propagating  on  a  current  may  experience  blocking,  if  the  current 
is  sufficiently  strong  and  has  a  component  opposing  the  waves.  If  blocking  occurs,  the  wave 
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energy  cannot  be  transported  against  the  current,  and  the  waves  are  “stopped”  (in  the  real  case 
this  often  implies  confused  wave  conditions  that  could  cause  navigational  hazards).  The 
criterion  for  determining  the  limit  for  wave  blocking  is  given  by  (Jonsson,  Skougaard,  and  Wang 
1970): 

Cgr  +  U  cos((5  -  a)  =  0  (13) 


If  Equation  13  is  fulfilled,  the  denominator  in  Equation  9  becomes  0  and  the  wave  rays  (along 
which  the  energy  is  conserved)  are  perpendicular  to  the  wave  orthogonals.  At  the  point  of 
blocking,  the  wavelength  attains  its  minimum  value,  which  may  be  estimated  from: 


./—tanh  kd  = 
V  L 


1  IIT 

1  —  n]]  L0 


where 


1 

n  =  — 
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sinh2 kd  J 
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(15) 


The  required  blocking  speed  associated  with  Equation  14  may  be  estimated  from  Equation  13, 
once  the  wavelength  L  at  blocking  has  been  determined  for  a  specific  La  and  d.  This  criterion 
may  be  written  in  non-dimensional  form  as: 


U  cos(8-a)  n  L  d  ,  ,  , 

- - -  =  - — - tanh  kd 

§Ta  2k  Vd  Lo 


Thus,  for  a  specific  ratio  d/L0  the  required  blocking  speed  might  be  determined  from 
Equations  14  and  16.  Figure  2  displays  the  non-dimensional  blocking  speed  as  a  function  of 
d/L0. 


Asymptotic  solutions  to  the  conditions  for  blocking  may  be  readily  obtained  for  shallow  and 
deep  water.  In  deep  water,  that  is,  kd  — »  °°,  Equation  14  yields: 

L-—L  (17) 

4  °  v  J 

and  Equation  16  produces: 

U  cos(S-a)  =~~Ca  (18) 
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Figure  2.  Non-dimensional  blocking  speed  as  a  function 
of  relative  water  depth 


where  C0  is  the  phase  speed  in  deep  water  neglecting  the  current  ( '=gTJ2n ).  In  shallow  water,  kd 
becomes  small;  omitting  terms  of  order  (kd)2  and  higher  results  in  the  following  expression  for 
the  wavelength  at  blocking: 


d_ 

L 


( 


9  d 
327i5  T 


\l/6 


(19) 


The  corresponding  current  speed  at  blocking  is  given  by: 


U  cos(S-a) 

gTa 


1  ( d  ' 
L 


1/2 


(20) 


This  relationship  is  in  fact  identical  to  U  cos(8-a )=(gd)112.  The  asymptotes  for  deep  and  shallow 
water  are  indicated  in  Figure  2. 

The  wave  equations  are  numerically  solved  by  a  finite-difference  formulation  similar  to  that  in 
NMLONG.  For  representing  waves  with  random  height,  a  Monte-Carlo  simulation  is  carried  out 
for  a  large  number  of  individual  waves  belonging  to  the  Rayleigh  probability  density  distribution 
in  the  offshore  so  that  statistically  stable  wave  quantities  are  obtained  in  averaging  the  results 
from  all  the  waves.  The  Rayleigh  distribution  will  not  necessarily  be  followed  as  the  waves 
propagate  onshore  and  break  or  are  modified  by  the  current-wave  interaction.  Because  the  wave 
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transformation  depends  on  the  current  field,  the  wave  calculations  are  updated  after  the  current 
has  been  computed  (described  next).  This  iteration  between  the  wave  and  the  current  field  is 
continued  until  convergence  is  achieved. 

LONGSHORE  CURRENT:  After  the  wave  transformation  calculations,  the  longshore  current 
is  computed  using  the  alongshore  momentum  equation  where  lateral  mixing,  bottom  friction,  and 
external  forcing  are  included: 


d  (  ,dV 

—  8  d - 

dx  I  dx 


ft 


by 


1  dS^ 

p  dx 


-R..-R, 


(21) 


where  V  is  the  longshore  current  velocity,  the  bottom  friction  stress  (discussed  later),  8  a 
lateral  mixing  coefficient  [(=A Hum,  where  H  is  the  wave  height,  for  random  waves  taken  to  be 
the  root-mean- square  (rms)  wave  height,  um  the  bottom  orbital  velocity,  and  A  an  empirical 
coefficient  typically  in  the  range  0. 2-0.5),  Sxy  the  radiation  stress  transported  onshore  and 
directed  alongshore,  and  Rw  and  R/c  forcing  associated  with  wind  and  an  external  (large-scale) 
current  (e.g.,  tide),  respectively.  The  velocity  V  constitutes  the  alongshore  component  of  U\  that 
is,  U=(V2+U2)U2,  where  Uc  is  the  mean  cross-shore  velocity  (thus,  the  angle  8  is  given  by  tan' 
\V/Uc)l 

The  forcing  associated  with  a  local  wind  is  determined  by: 


RW=CD^\W\W  sirup  (22) 

P 

where  Cd  is  a  drag  coefficient  given  by  the  expression  developed  by  the  WAMDI  Group  (1988), 
pa  the  air  density,  W  the  wind  speed,  and  cp  the  wind  direction  (W  and  cp  are  defined  in  the  same 
way  as  the  current;  see  Figure  1).  In  NMLong-CW,  it  is  possible  to  specify  an  external  current, 
assumed  to  be  associated  with  a  large-scale  circulation,  such  as  the  tide  or  a  regional  coastal 
current.  To  represent  this  current  in  the  model,  the  forcing  is  derived  from  a  term  introduced  as: 

Rlc=Cf\Ulc\Ulc  (23) 

where  cyis  the  bottom  friction  coefficient  appearing  in  fby  (typically  in  the  range  0.002-0.008  for 
field  conditions),  and  Uic  the  specified  longshore  component  of  the  external  current  (the  cross¬ 
shore  component  of  this  current  is  equal  to  Uc ).  The  bottom  friction  stress  fby  is  calculated  with 
Nishimura’s  square-wave  approximation  (Nishimura  1988)  as  described  by  Kraus  and  Larson 
(1991)  to  save  substantial  execution  time  yet  incorporate  the  non-linear  term. 

Finally,  the  radiation  stress  Sxy  is  calculated  from: 

Sxv=^-pgH2^sm2a  (24) 

to  Cr 
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MEAN  WATER  LEVEL:  The  mean  water  level  T|  (setup  and  setdown)  is  determined  from  the 
cross-shore  momentum  equation: 


dr i  dSxx  I  , 

—  = - —  -  CDpa  \W  W  cos  tp 

dx  dx 


(25) 


in  which  Sxx  is  the  radiation  stress  component  transported  and  directed  onshore,  given  by: 


(26) 


For  random  waves,  Sxx  and  Sxy  are  determined  as  averages  for  the  selected  number  of  waves  in 
the  Monte-Carlo  simulation  before  they  are  inserted  in  the  momentum  equations.  The  cross¬ 
shore  mean  current  does  not  enter  Equation  25  and  is  assumed  only  to  modify  the  wave 
transformation. 

EXAMPLES: 

Example  1:  Wave  transformation  at  an  inlet  entrance.  Smith  et  al.  (1998)  measured 
wave  breaking  on  a  current  at  an  idealized  inlet  in  the  laboratory.  A  1:50  scale  model  of  an  inlet 
was  constructed  in  a  46-m-wide  by  99-m-long  concrete  basin  with  0.6-m-high  walls.  The 
parallel  jetties  at  the  inlet  had  a  spacing  of  3.66  m  and  extended  5.5  m  offshore.  A  seaward 
flowing  (ebb)  current  Uc  was  generated  between  the  jetties  that  diffused  as  it  propagated 
offshore.  The  experimental  conditions  constituted  permutations  of  the  following  parameter 
values:  mean  spectral  wave  height  in  deep  water  Hmo  =3.7  and  5.5  cm,  spectral  peak  period  Tp  = 
0.7  and  1.4  sec,  wave  direction  perpendicular  to  the  jetties,  and  Uc  =  0,  12,  and  24  cm/sec.  Wave 
height  and  current  were  measured  at  several  gauges  around  the  inlet  with  the  main  objectives  to 
study  wave  breaking  and  determining  the  decay  in  wave  height  on  the  current. 

Here,  two  cases  are  discussed  to  illustrate  the  performance  of  NMLong-CW,  in  particular  for  the 
algorithm  developed  to  calculate  wave  breaking  on  a  current.  These  simulations  were  partly 
carried  out  to  validate  the  generalization  of  Equation  8  (calculation  of  the  energy  dissipation  on  a 
current).  It  was  found  that  the  standard  value  of  V=  0.4  overall  provided  reasonable  results, 
although  this  value  should  be  confirmed  by  further  simulations  against  other  data  sets.  The  cases 
discussed  here  encompassed  Case  5  (Hmo  =  4.1  cm,  Tp  =  1.4  sec,  Uc=  13  cm/sec)  and  Case  11 
(Hmo  =  3.9  cm,  Tp  =  0.83  sec,  Uc  =  24  cm/sec)  from  Smith  et  al.  (1998)  illustrating  the  results  for 
both  the  weaker  and  stronger  current  cases.  Wave  height  transformation  was  calculated  by 
Monte-Carlo  simulation  with  a  Rayleigh  distribution  specified  in  the  offshore.  The  actual 
simulated  time  series  of  waves  at  each  location  was  input  in  the  present  cases  to  compute  the 
significant  wave  height,  which  was  assumed  to  be  equal  to  Hmo.  Figures  3a  and  3b  display  the 
calculated  significant  wave  height  for  Cases  5  and  11,  respectively,  together  with  the  measured 
wave  heights.  Viewing  Figure  3  together  with  other  cases  not  shown  here,  the  agreement 
between  calculations  and  measurements  is  regarded  as  satisfactory. 
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(a) (b) 

Figure  3.  Comparison  between  calculated  and  measured  significant  wave  height  for 
Smith  etal.  (1998)  (a)  Case  5  and  (b)  Case  11 

Example  2:  Wave-generated  longshore  current  in  the  presence  of  a  large-scale 
current.  To  illustrate  the  capability  in  NMLong-CW  to  simulate  the  presence  of  a  large-scale 
current  and  the  wave-generated  longshore  current  in  the  nearshore,  a  hypothetical  example  is 
discussed.  An  equilibrium  profile  shape  in  accordance  with  Dean  (1977)  was  assumed  with  a 
shape  parameter  of  A  =  0.1  m1/3  corresponding  to  a  median  sand  grain  size  of  about  0.2  mm.  A 
root-mean- square  (rms)  wave  height  in  deep  water  of  Hrmw  =  2.0  m  with  a  mean  period  of  T 
=  8.0  sec  and  a  mean  incident  angle  a0  =  30  deg  was  specified  (waves  Rayleigh  distributed  in 
deep  water).  Also,  a  large-scale  current  was  specified  with  an  alongshore  component  growing 
exponentially  from  0  at  the  shoreline  to  0.5  m/sec  in  the  offshore  and  having  no  component 
across  shore  (8  =  90  deg).  Standard  (default)  values  were  specified  for  the  wave  and  longshore 
current  parameters:  y ),  =  0.78,  K  =  0.15,  T  =  0.4,  A  =  0.3,  and  c/=  0.003. 

Figure  4  illustrates  the  simulated  longshore  current  for  waves  and  large-scale  (L-S)  current 
together  as  well  as  for  waves  only.  Also,  the  cross-shore  distribution  of  the  input  L-S  current  is 
shown  in  the  figure.  In  the  absence  of  waves,  NMLong-CW  exactly  reproduces  the  input  L-S 
current.  However,  if  waves  are  present,  the  bottom  friction  stress  will  increase,  and  the 
simulated  current  will  typically  also  differ  from  the  L-S  current  outside  the  region  of  wave¬ 
generated  currents.  Because  the  influence  of  the  waves  disappears  in  deeper  water,  the  simulated 
current  will  approach  the  input  L-S  current  (if  there  is  no  wind-induced  current). 

Example  3:  Wave  transformation  at  an  inlet  in  the  presence  of  a  flood  and  ebb 
current.  This  example  involves  inlet  currents  and  waves  representative  of  hydrodynamic 
conditions  as  observed  at  Shinnecock  Inlet,  Long  Island,  NY.  For  this  example,  NMLong-CW  is 
operated  to  obtain  information  exclusively  on  cross-shore  processes.  It  is  cautioned  that  inlet 
entrances  are  complex  and  that  results  obtained  from  simplified  calculation  conditions  should  be 
interpreted  with  caution. 
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Figure  4.  Simulation  of  the  effect  of  a  large-scale  current 
on  a  wave-generated  nearshore  current 


The  example  concerns  wave  transformation  on  flood  and  ebb  currents.  An  inlet  channel  is 
simulated  with  a  water  depth  of  4  m  at  the  throat  and  linearly  sloping  offshore  to  a  depth  of  12  m 
(assumed  boundary  for  this  inlet).  Waves  are  assumed  to  travel  along  the  channel  (0  incident 
wave  angle),  and  the  tidal  current  (flood  or  ebb)  decreases  linearly  from  the  throat  to  the  offshore 
end  of  the  channel  where  it  was  set  to  0  (at  the  12-m  water  depth  corresponding  to  x  =  0  in 
Figure  5).  The  deepwater  rms  wave  height  was  1.0  m,  and  the  mean  wave  period  8  sec  (typical 
for  Shinnecock  Inlet).  Standard  parameter  values  were  employed  in  the  NMLong-CW 
simulations. 

Figure  5  displays  the  results  of  the  simulations  for  two  different  current  speeds  at  the  inlet  throat 
for  flood,  ebb,  and  no-current  cases.  In  the  no-current  case,  a  small  increase  in  wave  height  is 
observed  because  of  shoaling,  whereas  for  the  flood  current  the  waves  experience  a  reduction  in 
height  as  they  approach  the  inlet.  The  opposite  trend  occurs  if  the  waves  encounter  an  ebb  flow 
and  a  pronounced  wave  height  increase  might  take  place.  For  example,  in  the  case  of  an  ebb 
current  of  3.0  m/sec  at  the  inlet  throat,  wave  breaking  occurs  because  of  the  limitation  in  wave 
steepness,  causing  a  reduction  in  wave  height,  as  seen  in  Figure  5  (note  that  the  waves  propagate 
along  the  x-axis  to  allow  qualitative  comparison  with  Figure  6).  Also,  for  the  larger  ebb  current, 
wave  blocking  is  taking  place  before  the  waves  reach  the  inlet  throat.  Figure  6  is  a  photograph 
taken  from  the  east  jetty  at  Shinnecock  Inlet  and  shows  wave  breaking  and  blocking  on  an  ebb 
current.  Note  that  the  waves  cannot  penetrate  against  the  current,  with  turbulent  water  to  the  left 
(south)  and  calm  water  to  the  north,  inside  the  inlet. 
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Figure  5.  Calculated  wave  transformation  on  flood  and  ebb 
currents  at  an  inlet 


Limit  of  wave  penetration 


Figure  6.  Wave  breaking  and  blocking  by  an  ebb  current 
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By  this  example,  it  can  be  seen  that  the  wave  climate  in  an  inlet  channel  can  be  investigated  at 
reconnaissance  level  with  NMLong-CW  under  the  assumption  of  longshore  uniformity.  For 
example,  under  a  given  ebb  current  and  offshore  wave  height  and  period,  the  increase  in  wave 
height  and  in  wave  steepness,  defined  as  H/L ,  owing  to  the  presence  of  the  tidal  current  can  be 
calculated.  Steep  waves  pose  a  hazard  to  navigation  if  the  wavelength  approaches  that  of  the 
vessel  transiting  the  inlet. 

FUTURE  WORK:  Research  is  underway  in  the  CIRP  to  include  longshore  sediment  transport 
in  NMLong-WC.  The  model  and  interface  will  be  revised  to  calculate  longshore  transport  inside 
and  seaward  of  the  surf  zone  as  associated  with  a  current  produced  by  waves,  tide,  and  wind. 

ADDITIONAL  INFORMATION: 

This  CETN  was  written  by  Dr.  Magnus  Larson  of  the  Department  of  Water  Resources 
Engineering,  Lund  Institute  of  Technology,  University  of  Lund,  Sweden,  and  by  Dr.  Nicholas  C. 
Kraus  of  the  U.S.  Army  Engineer  Research  and  Development  Center  (ERDC),  Coastal  and 
Hydraulics  Laboratory.  The  work  performed  at  the  University  of  Lund  was  done  under  contract 
with  the  ERDC  through  the  U.S.  Army  Engineer  Research  and  Standardization  Group-United 
Kingdom.  Questions  about  this  CETN  can  be  addressed  to  Dr.  Nicholas  C.  Kraus  (601-634- 
2016,  Fax  601-634-3080,  e-mail:  KrausN@wes.army.mil). 

This  CETN  should  be  cited  as  follows: 

Larson,  M.,  and  Kraus,  N.C.,  (2000).  “Enhancements  of  the  numerical  model  of  the 
longshore  current  NMLONG  to  include  interaction  between  currents  and  waves 
(NMLong-CW),”  Coastal  Engineering  Technical  Note  CETN-IV-25,  U.S.  Army 
Engineer  Research  and  Development  Center,  Vicksburg,  MS, 
http://chl.wes.army.mil/library/publications/cetn/. 
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